PROCEEDINGS 

°F SCIENCE 



p meson decay width from 2+1 flavor lattice QCD 



o 

(N 

> 
O 

(N 



> 

m 
o 



N. Ishizuka*" for PACS-CS Collaboration, 

" Graduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba 305-8571, 
Japan 

^ Center for Computational Sciences, University of Tsukuba, Tsukuba 305-8577, Japan 



We perform a lattice QCD study of the p meson decay from the Nf = 2+1 full QCD configu- 
rations generated with a renormalization group improved gauge action and a non-perturbatively 
(9 (a) -improved Wilson fermion action. The resonance parameters, the effective p — > nn cou- 
pling constant and the resonance mass, are estimated from the P-wave scattering phase shift 
for the isospin 7=1 two-pion system. The finite size formulas are employed to calculate the 
phase shift from the energy on the lattice. Our calculations are carried out at two quark masses, 
= 410MeV (m^/mp = 0.46) and ;7% = 300 MeV (m,,/mp = 0.35), on a 32^ x 64 (La = 2.9 fm) 
lattice at the lattice spacing a = 0.091 fm. We compare our results at these two quark masses with 
those given in the previous works using Nf = 2 full QCD configurations and the experiment. 
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1. Introduction 



Study of the p meson decay is a significant step for understanding the dynamical aspects of 
hadron interactions with lattice QCD. In the early stage of studies toward this direction the transi- 
tion amplitude (?r7r|p) extracted from the time behavior of the correlation function {n{t)n{t)p{0)) 
was used to estimate the decay width, assuming that the hadron interaction is small |jl|, ^, Q]. 

A more realistic approach is a study from the P-wave scattering phase shift for the isospin 7=1 
two-pion system. The finite size formulas presented by Liischer in the center of mass frame [||] and 
extensions to non-zero total momentum frames ||^, ^ are employed for an estimation of the phase 
shift from an eigenvalue of the energy on the lattice. The first study of this approach was carried out 
by CP-PACS Collaboration using Nf = 2 full QCD configurations {m^ = 330 MeV, a = 0.2 Ifm, 
La = 2.5 fm) [^]. After this work ETMC Collaboration presented results with Nf = 2 configurations 
at several quark masses {m^ = 290,330MeV (La = 2.5 fm), = 420,480MeV (La = 1.9fm), 



a = 0.079 fm) [M 10]. Recently Lang et al. reported results of high statistical calculations on a 



single Nf = 2 gauge ensemble (m^r = 266MeV, a = 0.124fm, La = 1.98 fm) 

In the present work we extend these studies by employing Nf = 2 + 1 full QCD configurations 
and working on a larger lattice volume. Our calculations are carried out with the gauge config- 
urations previously generated by PACS-CS Collaboration with a renormalization group improved 
gauge action and a non-perturbatively O (a) -improved Wilson fermion action at j3 = 1 .9 on 32^ x 64 
lattice (a = 0.091 fm. La = 2.9fm) We choose two subsets of the PACS-CS configurations. 
One of them corresponds to the hopping parameters K"„^ = 0.13754 for the degenerate up and 
down quarks and = 0. 13640 for the strange quark, for which the pion mass takes niji = 410MeV 
(m:t/mp = 0.46). The other is at K^d = 0. 13770 and = 0.13640, corresponding to m^j = 300 MeV 
{mji/nip = 0.35). All calculations are carried out on the PACS-CS computer at Center for Compu- 
tational Sciences, University of Tsukuba. Details of our calculations are presented in Ref. 1 13]. 

We note that QCDSF Collaboration calculated the scattering phase shifts for the ground state 
in the center of mass frame at several quark masses, {rrij^ = 240 — 430 MeV) [Q. They estimated 
the resonance parameters from these results, assuming that the effective p — )■ nn coupling constant 
does not depend on the quark mass. BMW Collaboration presented their first preliminary results 
with Nf = 2 + 1 configurations (m^j = 200, 340 MeV, a = 0.116fm) at Lattice 2010 [|E 



2. Method 



In order to calculate the f -wave scattering phase shift for the isospin 7=1 two-pion system at 
various energies from a single full QCD ensemble, we consider three momentum frames, the center 
of mass frame (CMF), the non-zero momentum frames with total momentum P = (27r/L)(0,0, 1) 
(MFl) and P = (27r/L)(l, 1,0) (MF2), as carried out in the previous works by ETMC [g, |l^ and 
Lang et al. [pi]]. In the present work we calculate the scattering phase shifts for four irreducible 
representations : Tj" in the CMF, A,^ and in the MFl, and Bj" in the MF2. The finite size 
formulas for these representations are given in Refs. [Q, ^, ^]. 

For the Tj^ and the representation, we only calculate the scattering phase for the ground 
state in the present work. The energy of the ground state is much smaller than that of the excited 
state on our gauge configurations. Thus the energy of these states can be extracted by a single 
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exponential fit for the time correlation functions of the p meson. We use the local p meson operator 
for the sink and a smeared operator for the source as discussed later. 

For the A2 and the Bj" representation, we also calculate the scattering phase shift for the first 
excited state. In order to extract the energies of the lowest two state for these representations, we 



use the variational method [17] with a matrix of the time correlation function. 



^ , (oi [nKy{t){7L7Lm io) (oi {KKymits) io) , 

' (0|pt(0(7r7r)(O|0) (0|pt(Op(^.,)|0> 



for each representation. The energies are extracted from two eigenvalues A„(f) {n= 1,2) of the ma- 
trix M{t) = G{t) {tR) with some reference time Ir, assuming that the lower two states dominate 
the correlation function. 

In (lA) the operator p{t) is given by p{t) = Y^j=i Pj ■ P7(p,0/|Pl' where Pj{p,t) is the local 



operator for the neutral p meson at the time slice t with the momentum p. The momentum takes p = 
(27r/L)(0,0, 1) for the and p = (27r/L)(l, 1,0) for the Bj" representation. Hereafter we assume 
that the momentum p takes one of these two values depending on the representation. {mi){t) is an 
operator for the two pions with the momentum and p, which is defined by 

(7r7r)(0 = ^(7r+(O,fi)7r-(p,0-7r-(O,?i)7r+(p,0) xe^-C'-') , (2.2) 

where 7r^(p,f) is the local pion operator with the momentum p at the time slice t. The time slice of 
the pion with the zero momentum is fixed at t\ » t, and the time slice of the other pion t runs over 



the whole time extent. An exponential time factor in (2.2) is introduced so that the operator has the 
same time behavior as that of the usual Heisenberg operator, i.e., (0| (7r7r)^(f) = (0| (7171)^(0) e^^', 
with the Hamiltonian H. 



Two operators {KK){ts) and p(f.v) are used for the sources in (2.1 ), which are given by 



{KK){t,) = -^(71+ (0,?,)7r"(p,f.) -71^(0, f,)7r+(p,f,)) , (2.3) 

P^^x) = E 4^ {u{i,h)YpU{i,t,) -D(z,f,)7p£>(z,f,))e'P-^ , (2.4) 
V 2 ^ ^ 

where Jp = Ys'j=\_Pj " T/'/IpI- The operator Q{z,ts) (Q = U,D) is a smeared operator for the up or 
the down quark given by Q{z,ts) = L^^(x,fv) -^^dx — z|), where q{x,ts) (q = u,d) is the up or the 
down quark at the position x and the time tg. We adopt the same smearing function ^I'dx]) as in 
Ref. [p^]. This operator is used after fixing gauge configurations to the Coulomb gauge. In ( ^^ 
a summation over z is taken to reduce a statistical error and F = { z | z = (L/2) • (ni , ?i2> '^3) 1 '^y = 
or 1 , A^r = 8 } is chosen in the present work. The smeared operator ( [2.4[ ) is also used to extract 
the energy of the ground state for the and the representation, setting the momentum p = 
and p = (27r/L)(0,0, 1), respectively. 

The periodic boundary conditions are imposed for both spatial and temporal directions in 
configuration generations. We impose the Dirichlet boundary condition for the temporal direction 
in calculations of the quark propagators, to avoid the unwanted thermal contributions produced by 
propagating two pions in opposite directions in a time. For both quark masses, we set the source 
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Figure 1: Effective masses for six states considered in the present work at irijc = 410MeV (left panel) and 
;7% = 300 MeV (right panel). 



operators at ts = \2 to avoid effects from the temporal boundary, and the zero momentum pion in 
the sink operator {KK){t) in ( |Z2| ) at ti = 42. The total number of configurations at niji = 410MeV 
is 440 and that nin = 300 MeV is 400. We calculate the quark contractions of G{t) in ( |01 ) by the 
source method and the stochastic noise method as in the previous work by CP-PACS [||]. Details 
of the method of the calculations are explained in Ref. |[T3|]. 



3. Results 



The effective masses of the time correlation functions for six states considered in the present 
work are plotted in Fig. |l|, where we choose tn = 23 as the reference time of the variational method. 
We can find plateaus in the time region t >23. We extract the energy by a single exponential fit 
for the time correlation functions. In Fig. [T] the results of the fitting with one standard deviation 
error band are also expressed by solid lines. The dotted line for the and representation in 
the figure indicates the energy of the two free pions for each representation. 
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Figure!: {k^ /tan5{k))/y^ as a function of square of the invariant mass s at m^i = 410MeV (left panel) 
and mj[ = 300 MeV (right panel). 



Converting the energies for each representation to the invariant masses ^/s and substituting 
them into the finite size formulas, we obtain the scattering phase shifts plotted in Fig. |2[ In the figure 
we show {k^ / tan 5{k))/ y/s as a function of square of the invariant mass s, where k = ys/A^ 



-m'l 



is the scattering momentum. The finite size formulas for the and the Bj^ representation are 
plotted by dotted lines. 

In order to extract the resonance parameters from our results of the scattering phase shift, we 
parametrize the phase shift with the effective p — tttt coupling constant gpnn and the the resonance 
mass YHp by 

-477tM=4^K--). (3.1) 

This parametrization has been widely used in the previous works of the p meson decay. By chi- 



square fitting of the scattering phase shifts with the fit function (3. 1 ), we obtain. 



,;r;r= 5.52 ±0.40 
,;,;j= 5.98 ±0.56 



mp = 892.8 ±5.5 ±13MeV 
mp = 863 ±23 ±12MeV 



atm;r = 410MeV , 
atm;r = 300MeV , 



(3.2) 
(3.3) 



where the second error of nip is the systematic uncertainty for the determination of the lattice 
spacing. In Fig. § we draw the fitting curves by solid red lines. 



From ( |3.2[ ) and ( |3.3| ) we find that the gp%n at the two quark masses are consistent within the 
statistical error and also with the experiment gpnn = 5.874 ±0.014 given from the experimental 
results of the decay width Fp = 146.2 ± 0.7 MeV [18|. This suggests a weak quark mass depen- 



dence of the coupling constant. But our calculations are carried out only at the two quark masses, 
so calculations at more quark masses are necessary to obtain a definite conclusion for the quark 
mass dependence. We leave this issue to studies in the future. 

In Fig. |3| we compare our results (PACS-CS) obtained in 2 ± 1 flavor QCD with those by 
ETMC [^, and Lang et al. [ pi] ] in 2 flavor QCD. The left panel shows the effective coupling 
constant gp^K and the right panel displays the resonance mass lUp as a function of m\. A good 
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Figure 3: Comparison of our results (PACS-CS) obtained in 2 + 1 flavor QCD with those by ETMC and 
Lang et al. in 2 flavor QCD. Left panel shows the effective coupling constant gpjiK and right is the resonance 
mass mn- 



agreement between our result and ETMC is observed for gpjui. The result for the coupling constant 
by Lang et al. takes a sUghtly smaller value, but it is almost consistent with other works. 

We see, however, large discrepancy for the resonance mass mp in the right panel of Fig. ||. 
The three groups worked at a single lattice spacing, therefor a possible reason of the discrepancy is 
the discretization error due to the finite lattice spacing. We can also consider several other reasons, 
the reliability of the determination of the lattice spacing, the dynamical strange quark effect, the 
isospin breaking effect and so on, but a definite conclusion can not be given here. 



4. Summary 

We have reported on a calculation of the P-wave scattering phase shift for the isospin 1 = 1 
two-pion system and an estimation of the resonance parameters from the Nf = 2-'rl full QCD 
configurations with a large lattice volume. The calculations are carried out at two quark masses, 
which correspond to rrijc = 410MeV and 300MeV. Our results of the effective coupling constant 
gpTiK at the two quark masses are consistent within the statistical error and also with the experiment. 
This suggests a weak quark mass dependence of the coupling constant. We find a discrepancy for 
the resonance mass nip among three lattice studies. Calculations near or on the physical point 
closer to the continuum limit are necessary for a precise determination of the resonance mass from 
lattice QCD. We leave this issue to studies in the future. 
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